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The character of critical behavior in physical systems depends on the range of interactions. In the 
limit of infinite range of the interactions, systems will exhibit mean-held critical behavior, i.e., critical 
behavior not affected by fluctuations of the order parameter. If the interaction range is finite, the 
critical behavior asymptotically close to the critical point is determined by fluctuations and the actual 
critical behavior depends on the particular universality class. A variety of systems, including fluids and 
anisotropic ferromagnets, belongs to the three-dimensional Ising universality class. Recent numerical 
studies of Ising models with different interaction ranges have revealed a spectacular crossover between 
the asymptotic fluctuation-induced critical behavior and mean-field-type critical behavior. In this 
work, we compare these numerical results with a crossover Landau model based on renormalization- 
group matching. For this purpose we consider an application of the crossover Landau model to the 
three-dimensional Ising model without fitting to any adjustable parameters. The crossover behavior 
of the critical susceptibility and of the order parameter is analyzed over a broad range (ten orders) 
of the scaled distance to the critical temperature. The dependence of the coupling constant on the 
interaction range, governing the crossover critical behavior, is discussed 



I. INTRODUCTION 

In systems with a critical-point phase transition the 
asymptotic critical behavior of thermodynamic properties 
can be characterized by scaling laws with universal critical 
exponents and universal scaling functions jl| . In this issue 
of the Journal of Statistical Physics, dedicated to M. E. 
Fisher, it is interesting to recall some of the discussions at 
a conference on critical phenomena held in Washington, 
DC in 1965 At that conference, Fisher reviewed the 
definitions of the proposed asymptotic critical power laws, 
both for the thermodynamic properties B a nd for the cor- 
relation function of the order parameter Commenting 
on the presentations of Fisher, Debye || made the follow- 
ing remark: "I would like that the theoretical people tell 
me when I am so and so far away from the critical point, 
then my curve should look so and so," to which Fisher re- 
sponded (|]: "What Dr. Debye said is correct, namely that 
when one is talking about some deviations one must at- 
tach to it a magnitude and a range of temperatures where 
it is to be observed, which can be different for one phe- 
nomenon compared to another. When these ranges differ 
from one system to another for the same phenomenon, we 
have a real challenge." It is the purpose of the present pa- 
per to further elucidate the issue of nonasymptotic critical 
behavior. 

The asymptotic fluctuation-induced critical scaling laws 
are valid sufficiently close to the critical point, where the 
correlation length £ of the long-range critical fluctuations 
will be much larger than the range of the intcrmolccular 
interaction. Away from the critical point, the correlation 



length £ decreases and sufficiently far away from the crit- 
ical point £ will become of the order of the intermolecular 
interaction range and the effects of any critical fluctuations 
become negligible. 

The actual critical behavior will depend strongly on the 
range of the microscopic interactions. As the interaction 
range becomes larger and larger, the critical fluctuations 
become more and more suppressed and mean-field the- 
ory will become applicable up to temperatures closer and 
closer to the critical temperature. In the limit of infinite 
interaction one recovers classical, i.e., mean- field asymp- 
totic critical behavior. Hence, as the interaction range 
varies, the system will exhibit a crossover from fluctuation- 
induced critical behavior to mean-field critical behavior. 

Recently, Luijten and coworkers have reported numer- 
ical studies of the susceptibility and the order parame- 
ter of both two-dimensional [0-^) and three-dimensional 
[[Hi -12 Ising models with different interaction ranges, 
showing crossover from fluctuation-induced to mean-field 
critical behavior. On the other hand, during the past two 
decades theoretical equations for dealing with crossover 
from fluctuation-induced to mean-field critical behavior 
have been developed by several groups of investigators 
[|l3 29 1 . A comparison between some of these theoretical 
approaches ]l6|- f2"l| , p4[ has been presented by Anisimov et 
al. [Q. Here we consider a crossover Landau model that 
has been used to represent experimental thermodynamic- 
property data in the critical region of fluids, which are 
expected to belong to the three-dimensional Ising univer- 
sality class @U!oU3l). Thc 

crossover Landau model 
contains some system-dependent parameters related to 
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the coefficients in a Landau-Ginzburg Hamiltonian B2 ], 
Since these system-dependent coefficients are not known 
a priori for fluids, they are used as adjustable parame- 
ters in fits to experimental data. However, in the case 
of the fsing models studied numerically by Luijten and 
Binder we have a priori information concerning 

the system-dependent constants. Hence, it has become 
possible to evaluate the crossover Landau model for fsing 
models without using adjustable parameters and to make 
a comparison with the numerical data obtained by Luijten 
and Binder for three-dimensional Ising models. 



II. CROSSOVER LANDAU MODEL 

A procedure for constructing an expression for the 
Helmholtz-energy density for fluids in the critical region, 
that incorporates crossover from fsing-like critical behav- 
ior to classical behavior was originally proposed by Chen 
and coworkers 19-2^]. The procedure, based on a theo- 
retical analysis of Nicoll and coworkers [p~3| |l5|] , represents 
an approximate solution of the nonlinear renormalization- 
group equations and involves a transformation to a clas- 
sical (Landau) expansion. Two crossover Landau mod- 
els, a two-term and a six-term crossover Landau model 
have been considered depending on the number of terms 
retained in the classical Landau expansion as further re- 
viewed by Anisimov et al. p3| . Explicit expressions for the 
amplitudes of the critical power laws incorporated in the 
crossover Landau models have been evaluated by Tang et 
al. p2] and, more recently, by Agayan et al. (3^]. Here we 
consider the two-term crossover model which deals with 
the crossover from asymptotic Ising-like critical behavior 
to asymptotic mean-field critical behavior p2j. 

In the crossover Landau model one starts with an 
asymptotic Landau expansion of the critical part AA of 
the classical local Helmholtz free-energy density A |22] : 



AA = ±a Q T<j) 2 + ±u <j) 4 + 5C (V(/>) 2 



(2.1) 



where r = (T — T c )/T is the reduced difference between 
the actual temperature T and the critical temperature T c , 
4> = 4>{y) is a spatially fluctuating order parameter, while 
a , Wo aim c o are system-dependent coefficients. For the 
Ising model the order parameter <f> corresponds to the mag- 
netization density. It is convenient to transform the vari- 
ables and coefficients, such that p3[ 



M 
a 



c P <t>, 
c 2 „c t , 



t = ctr, V = («d /3 A)V, (2.2) 
u*uAc% c = c 2 p v 2 /3 /7r 2 , (2.3) 
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where vq is the volume of an elementary lattice cell, 
w* ~ 0.472 is the universal coupling constant at the Ising 
fixed point |^2|,|4), u is the scaled tp 4 coupling constant, A 
is a dimensionless cut-off wave number, while Cj and c p are 
amplitudes of the rescaled temperature scaling field and of 



the rescaled magnetization density, respectively. Expres- 
sion (2.1) for A A can then be written as 

i(VAf) 2 . 



AA = \tM 2 



±u*uAM 4 



(2.4) 



Upon applying a transformation derived from renorma- 
lization-group matching [ p0| , pl| one obtains a renormalized 
critical part of the Helmholtz-energy density: 

WW A , r A^o~ . 1 .9,^ /r . r \ 

(2.5) 



AA S = \tM 2 TV 



4! 



-M 4 V 2 U 



where, in practice, the following approximate expressions 
for the rescaling functions T, T>, U, and K, are used: 



J = y(2^-l)/A B 

U = Y U ' A % 



and 



IC = 



auA 



1). 



(2.6) 
(2.7) 
(2.8) 

(2.9) 



Here v ~ 0.630, r\ ~ 0.033, and a ~ 0.110 are the criti- 
cal exponents associated with the correlation length, the 
spatial dependence of the order-parameter at the critical 
point and the heat capacity, such that dv = 2 — a, where 
d = 3 is the dimensionality [ |35|-^7 |, while A s ~ 0.51 is a 
correction-to-scaling exponent |35[ 58) . In the asymptotic 
crossover theory, the crossover function Y is to be evalu- 
ated from the equation p2|] 



1 - (1 - u)Y 



(2.10) 



where the parameter n plays the role of the distance from 
the critical point and it is given by 



= tT- 



"AM 2 VU. 



(2.11) 



In zero ordering field, M — 0, the isothermal suscepti- 



bilityxtx^M^A/^ 2 ), 
has an expansion of the form 

x = r+t-^i + r+t A 



(d 2 A/dM 2 ) tC - 2 \ above T c 



+ •••), 



(2.12) 



and along the line of spontaneous magnetization (coexis- 
tence curve) M = M cxc , below T c (t < 0): 

rrKI As 



x 



r |ir(i 



+ •••). 



(2.13) 



The order parameter CXC along the "coexistence curve" 
below the critical temperature (r < 0) has the expansion 

0cxc = c-Hl^ = ±Bv\t\ f \l + fliltl^ +•••)• (2.14) 

In the above equations, 7 = (2 — 77)1/ ~ 1.239 and 
(3 = (dv — 7)/2 ~ 0.326 are again universal critical ex- 
ponents |bB-BJ, while Tq, Tf, B and B\ are system- 
dependent amplitudes. These system-dependent ampli- 
tudes are related to wA and the coefficients in the trans- 
formations (2.2) and (2.3) and they are given by the fol- 
lowing equations p2| , ^3|j39| ] . For the leading amplitudes 
one has 
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2(/3-l/2) 



So = b I g 1 (ao/^o) 1/2 , 
and for the correction-to-scaling amplitudes 

2A S 



rf-rf(£] (i-<0, 



2A S 



(1-u). 



(2.15) 
(2.16) 
(2.17) 

(2.18) 
(2.19) 



The coefficients g , b, g^, and b\ in the above equation 
have been evaluated by Tang et al. |3^] .[] They depend on 
the critical exponents and the coupling constant u* and are 
therefore again universal constants with the values [(33"|]39| : 



and 



g + ~ 0.871, g~ ~ 0.174, 6-2.05, 



gt ~ 0.610, bi ~ 0.531. 



(2.20) 



(2.21) 



The theoretical expression for the amplitude Tj" in Eq. 
(2.13) has not yet been evaluated. The values implied by 
the crossover Landau model for the corresponding univer- 
sal amplitude ratios are 



r+/r - = 5.0, Bi/r+ = o.87, 



(2.22) 



to be compared with the values Tq /Lq = 4.95 ± 0.15 and 
Bi/Tf = 0.90 ±0.21 |||. A phenomenological reformula- 
tion of the two-term Landau model in terms of a paramet- 
ric representation yielding slightly more accurate values 
for some universal amplitude ratios has been presented by 
Agayan et al. p3[ . 

Wc define a crossover temperature t x as 



or, equivalently, 



(t2A) s 



(l-u)- 1/A % (2.23) 



(g+)- l / A °(uA) 2 (l-u) 



-l/A. 



(2.24) 



In the infinite-cutoff approximation (A — > oo, u — > 0, while 
uA remains finite), r x becomes the so-called Ginzburg 



*In Ref. |S2j Tang et al. actually considered a variety of ap- 
proximants for the rescaling functions T, T>, U, and /C. The 
crossover Landau model considered here corresponds to the one 
designated by Tang et al. as crossover model II. 



number [£2|;[LiJ . In general, the crossover behavior is gov- 
erned by two system-dependent parameters, namely, the 
cutoff A and the coupling constant u. However, for the 
Ising model A is a constant equal to unity independent of 
the range of interaction and only the coupling constant u 
governs the crossover behavior. One can also show that in 
the asymptotic crossover Landau model the appropriately 
scaled susceptibilities, ^ = X ±T Z/^o > an d the scaled or- 
der parameter, <j> = (frr^ /Bo become universal functions 
of a single argument r/r x ~t/t x |22[ . 

After some algebra, one can obtain the following uni- 
versal expression for the crossover behavior of the scaled 
susceptibility x + above T c in zero ordering field as a func- 
tion of t/t x : 



(9l 



y(i-7)/ 



As 




(2.25) 



where the rescaled crossover function Y, can be obtained 
by solving 



l-Y = (g+T/T x )- 1 / 2 Y<- 1 -^/ 2A °. 



(2.26) 



The scaled susceptibility x~ an d the scaled order param- 
eter <j) below T c are to be obtained by imposing the con- 
dition {3A/dM)t — for the phase boundary and solving 
Eqs. (2.5)-(2.9) numerically. 

We note that the crossover Landau model as specified 
above is an asymptotic version of the crossover Landau 
model @,^l[ valid for A/k^> 1. This asymptotic version 
appears to be adequate for the temperature ranges covered 
by the numerical studies of Luijten and Binder. 

The crossover Landau model has been derived from 
the renormalization-group theory of critical phenomena 
by using the technique of renormalization-group matching 
[ p^lpfjjil]! . In this approach one considers the relation be- 
tween AA at a given cutoff A and AA at another cutoff 
Ae~ l , where I is a variable. One then tries to determine 
a special "matching-point" value I = I* for which AA re- 
duces to the known classical expression in the absence of 
fluctuations. One limitation is that the value of I* that 
satisfies the matching condition for AA will not exactly 
satisfy the matching condition for its derivatives like X- 
Secondly, for practical reasons I* has been selected to sat- 
isfy the matching condition M = |2(],|2l|]. A comparison 
with the numerical data for the three-dimensional Ising 
model will enable us to make an assessment of the quality 
of these approximations. 
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III. APPLICATION TO THE 
THREE-DIMENSIONAL ISING MODEL 

In the Monte Carlo simulations a three-dimensional 
equivalent-neighbor model was adopted for which the 
Hamiltonian TL is given by JO] 



H/k B T = -J2 K ( r * 

(ij) 



rj)siS 



j - 



(3.1) 



where fee is Boltzmann's constant, Si = ±1, the sum runs 
over all spin pairs, and the spin-spin coupling is defined 
as K (r) = J > for |r| < R m and K{r) = for |r| > 
R m . However, in a renormalization-group analysis, various 
critical properties have been obtained as functions of an 
effective interaction range R defined by 



3+1 



with |r$ — Tj | < R Tl 



(3.2) 



where Kij stands for if(r,— Tj) while q is the coordination 
number. In the limit of an infinite interaction range R — > 
oo, R is related to R m by R 2 — 3i?^,/5, while R = R m for 
i? m = 1. The dependence of the effective interaction range 
R on R m is shown in Fig. 1. Also shown is an approximate 
relation El 



R 2 = irki 



lR- 2 ) 



(3.3) 



encompassing the known behavior in the limits R — » oo 
and R = 1. 

In order to apply the theoretical crossover model to 
the three-dimensional Ising model we need to specify the 
system-dependent constants. From the mean-field theory 
of the Ising model one obtains for the coefficients in the 
classical expansion (2.1) W3[ 



a = 1, 



uq = 2. 



(3.4) 



The physical cutoff wave number for the Ising model is 
7r/a, so that 



A = 1. 



(3.5) 



In a previous paper |42| ] , A of the three-dimensional Ising 
model was taken to be tt. However, we have subsequently 
realized that the spatial derivatives in the gradient term 
are already scaled by ix j a. 

A renormalization-group analysis of Eq. (3.1) has shown 
that u oc R~ A for large R |7|JTl]] . It then follows from Eq. 
(2.3) that ct oc R~ 2 and c p oc R and we write 



Up 

R 4 ' 



Ct 



cm 
i? 2 ' 



-1/2 



= Rc 



-1/2 



10 



(3.6) 



In a previous attempt to compare the crossover model 
with the numerical results for the three-dimensional Ising 



model, uq and c t o were treated as adjustable constants 
[ fT2| . In the present paper we shall adopt a priori esti- 
mates for uq and c t o- It has been suggested by Luijten 
[ fj"i"| that Mo j and consequently c t o, may actually exhibit a 
remaining weak dependence on R for small R. From Eqs. 
(2.3), (3.4) and (3.5) we note that the coefficients u (R) 
and Cto(R) are related by 



c 2 t0 (R) = lu*u (R). 



(3.7) 



The values of the coefficients uq and c t o can readily be 
determined for large R. For this purpose we note that the 
amplitude £g~ of the power law £ = £q~ t~ x / 2 for the clas- 
sical correlation length £ is given by £q — (co/ao) 1 / 2 p2| . 
Starting from an expression given by Fisher and Burford 
[ ffH and extending it to arbitrary interaction ranges, one 
finds that £q~ becomes equal to R/ V2d, so that Cq becomes 
equal to R 2 /2d for large R. Using Eqs. (2.3), (3.4) and 
(3.6) we thus obtain for large R: 



c t o(oo) = 6/tt 2 ~ 0.608, 



wo(oo) = — q (oo) ~ 1.566. 

(3.8) 



In the present paper we shall use two approximations 
for the coupling constant u of the three-dimensional Ising 
model. One estimate is obtained by assuming that u re- 
mains proportional to i?~ 4 for all values of R: 



U ((X)) 

R 4 



(3.9) 



Using Eq. (3.9) in conjunction with Eqs. (3.6)-(3.8) we 
can then predict from the crossover Landau model the am- 
plitude Tf which is related to the crossover temperature 
Tx through Eq. (2.23). The resulting values for as a 
function of R 3 are represented by the solid curve in Fig. 
2. We note that, since A s ~ 1/2, will approximately 
vary as R 3 for large R. 

A second estimate for Tj is obtained from numerical 
values calculated by Luijten for the amplitude b\ (R) of the 
correction-to-scaling contribution to the finite size effects 
[ jlj| . For large values of R where u approaches zero, b\{R) 
becomes proportional to (R)R ai ^ A " so that we may 
write 

T+(i?) = 6 (i?)&i(i?)i?~ Q7/As , (3-10) 

where bo(R) is a coefficient with a finite limiting value 

b (oo) = -2.6 ±0.5. (3.11) 

Values deduced for Tj from the numerical data obtained 
by Luijten for &i(i?), when the coefficient t>o(°o) m Eq. 
(3.10) is approximated by 6o(oo) for all R, are indicated 
by the circles in Fig. 2. These values can be reproduced 
from Eq. (2.23) for the crossover Landau model if we 
adopt the following approximation for u: 



4 



^o(oo) 
R 4 



(l + 0.9iT 



(3.12) 



as shown in the dotted curve in Fig. 2. 

Equations (3.9) and (3.12) are the two approximations 
that will be considered in applying the crossover Landau 
model to the Ising lattice. As we shall see, the predicted 
critical crossover behavior of the susceptibility is not very 
sensitive to these different choices for the coupling con- 
stant u for most values of R. 



IV. RESULTS OF ANALYSIS 

In Fig. 3 we present a plot of the scaled susceptibility 
X = X +T x/^o above T c as a function of t/t x . The solid 
curve represents the values predicted from Eq. (2.26) for 
the crossover Landau model with u — uq(oo)R~ in accor- 
dance with Eq. (3.9). The symbols indicate the numerical 
values obtained by Luijten and Binder (n]|llj] . A sensitive 
test of the crossover behavior is obtained by considering 
an effective susceptibility exponent 7^ defined by Q: 



its 



-dlogx/dlogr (r > 0). 



(4.1) 



Figure 4 shows the effective exponent 7^ as a function of 
r/r x . The symbols indicate values deduced by numerical 
differentiation of the susceptibilities obtained by Luijten 
]Tl| ] . The solid curve represents the effective exponent val- 
ues predicted from the crossover model. Note that this ex- 
ponent is related to the third derivative of the Hclmholtz- 
energy density AA. In Fig. 5 we show again 7^ as a func- 
tion of t/t x but now with u = -0 (oo)i?~ 4 (l + 0.9i? -2 ) in 
accordance with Eq. (3.12). We remark that the theoreti- 
cal prediction has a possible inaccuracy at the percentage 
level, since the values adopted for the various universal 
constants have a finite accuracy. In addition, we have 
only limited theoretical information for the values of u 
for small R. Based on the currently available theoreti- 
cal information we conclude from Figs. 4 and 5 that the 
crossover Landau model yields a satisfactory representa- 
tion of the crossover behavior of the susceptibility above 
T c for the three-dimensional Ising model with varying in- 
teraction ranges. 

As mentioned in Section 2, the crossover Landau model 
has been derived from the renormalization-group theory 
of critical phenomena by renormalization-group match- 
ing, but the match-point condition actually adopted corre- 
sponds to M = Q . This is correct for the susceptibility 
above T c at M = 0, but represents an approximation be- 
low T c where M cxc 7^ 0. Hence, we expect the crossover 
Landau model to be less accurate for a description of the 
critical crossover behavior of the susceptibility and of the 
order parameter below T c . 

In Fig. 6 we show the scaled order parameter <j) — 
4>t~^ /Bo and in Fig. 7 the scaled susceptibility xT = 



X~ T Z/^o as a function of — t/t x . The solid curves repre- 
sent again the values predicted from the crossover Landau 
model and the symbols indicate again the numerical data 
obtained by Luijten and Binder ]l0| j. In both cases the 
asymptotic Eq. (3.9) was adopted for u. While the ac- 
tual differences with the numerical values for </> and x~ 
are small and probably within the accuracy these proper- 
ties can be measured experimentally for real systems, the 
deviations are systematic. This becomes evident when one 
takes the logarithmic derivative so as to obtain the effec- 
tive order-parameter exponent 



P eS = dlog^/dlog |r|, (r < 0) 



(4.2) 



shown in Fig. 8 and the effective susceptibility exponent 
7 c ff 



7, 



efi 



-dlogx/dlog|r|, (t<0) 



(4.3) 



shown in Fig. 9. The effective exponents predicted by the 
crossover Landau model as a function of r/r x have the 
correct shape but the location of the crossover temperature 
r = t x is shifted from that implied by the numerical data. 
This shift amounts to about 1/2 of a decade for (3 e g and 
2/3 of a decade for 7^. The shift is larger for 7^ than for 
/3 e ff, since 7^ is related to the third derivative and /3 e ff is 
related to the second derivative of the crossover function 
deduced for AA by renormalization-group matching. 

Crossover Landau models with improved approximants 
for the rescaling functions T, V, U, and K, have been con- 
sidered by Tang et al. , but they never have been used 
for practical use because of increased complexity. 

The solid curves in Figs. 4, 5, 8 and 9 actually represent 
the effective exponents for R m > 1 where u < 1. For 
R m = R = 1 the leading amplitudes and Bq have 
been determined by Liu and Fisher [j35| from a careful 
analysis of series expansions. The value of = 1.095 
corresponding to 7 = 1.239 implies u = 1.17 as follows 
from Eq. (2.15). For u > 1, the correction-to-scaling 
amplitudes and B\ become negative in agreement with 
an earlier observation of Liu and Fisher . A comparison 
of the amplitudes calculated from the crossover model for 
u — 1.17 with the amplitudes deduced by Liu and Fisher 
| f?5[ from series expansions is presented in Table I. For 
u > 1, the asymptotic value of the effective susceptibility 
exponent 7^ is approached from above and that of the 
effective order-parameter exponent j3 e g from below |33],[l7| . 
The dependence of the effective exponents on r/r x with 
t x = |r^| _1 / As for the case R = 1, corresponding to 
u = 1.17, is represented by the dashed curves in Figs. 
4, 5, 8 and 9 in the temperature range covered by the 
numerical simulations. 

From Fig. 9 we note that the numerical data as well as 
the crossover Landau model show that 7^ passes through 
a minimum as a function of r/r x before reaching its mean- 
field value unity, (see Ref. 48). An interesting question, 
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raised earlier by Fisher [jl9), is whether also 7^ can as- 
sume values smaller than unity. Within the accuracy of 
our analysis, we do not find any evidence that the effec- 
tive susceptibility exponent above T c can assume values 
smaller than unity, as was already discussed in Ref. 10. 

V. DISCUSSION 

In the present paper we have applied a crossover Landau 
model, previously derived for the representation of ther- 
modynamic properties of fluids near the critical point to 
predict the critical crossover behavior of the susceptibil- 
ity above T c and of the order parameter and susceptibil- 
ity below T c for the three-dimensional Ising model with 
various interaction ranges. A comparison with numeri- 
cal data, recently obtained by Luijten and Binder [ fu^]i"l| , 
demonstrates that the crossover Landau model yields a 
good representation of the crossover behavior of the sus- 
ceptibility above T c . Below T c there are some systematic 
differences from the numerical data which are small for <j> 
and x~ themselves, but which become more pronounced 
when one considers the corresponding effective exponents. 
In applications of the crossover model to real fluids, resid- 
ual deviations are partially accounted for by using u and 
A, as well as the coefficients ao and uq in the classical ex- 
pansion (2.1) of the Helmholtz-energy density, or, equiva- 
lently, c* and c p (cf. Eq. (2.3)) as adjustable coefficients 

The analysis presented in this paper is based on an 
asymptotic version of the crossover Landau model which 
deals with the crossover from fluctuation-induced critical 
behavior to asymptotic mean-field critical behavior. For 
some other applications one needs to consider a more gen- 
eral version of the crossover Landau model that is ob- 
tained when the ratio A/re in Eq. (2.10) is replaced by 
(1 + A 2 /k 2 ) 1/2 fH). This generalization was originally in- 
troduced to account for crossover to nonasymptotic mean- 
field behavior in simple fluids f2l}| . Subsequently, it has 
become evident that this generalization becomes essen- 
tial in the case of complex fluids, when the crossover to 
mean-field behavior does not result from long-range inter- 
actions, but from a decrease of the effective cutoff wave 
number A. In that case A" 1 plays the role of a correla- 
tion length associated with an additional order parame- 
ter. Then the mean-field limit corresponds to mean- field 
tricriticality, when both £ and A -1 diverge |]50| , ^l| . 
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TABLE I. The critical amplitudes of the susceptibilities 
above and below T c and of the order parameter from the 
crossover two-term Landau model (CLM) using u — 1.17, and 
from series expansions [35] for R — 1. 
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FIG. 2. Amplitude of the correction-to-scaling contri- 
bution in the expansion (2.12) for the susceptibility. The 
solid curve represent values predicted from the crossover Lan- 
dau model with u = uo(oo)R~ 4 and the dotted curve with 
u — ito(oo)7? -4 (l + 0.9J? -2 ). The circles indicate numerical 
estimates deduced from finite-size effects obtained by Luijten 
[11]. The dashed lines represents the asymptotic linear depen- 
dence of T+ on R 3 . 



FIG. 1. Square of the effective interaction range R as a func- 
tion of R^- The symbols indicate values calculated by Lui- 
jten [11]. The solid curve represents the approximant given by 
Eq. (3.3). The dotted line represents the asymptotic relation 
R 2 = 3i?m/5 for large R. 




FIG. 3. Scaled susceptibility \ = X T Z^o as a function of 
r/r x . The symbols indicate numerical values from the Monte 
Carlo simulations [10,11]. The curve represents values calcu- 
lated from the crossover Landau model with u — Mo(oo)i? -4 . 
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FIG. 4. Effective susceptibility exponent 7^ as a function of 
t/tx- The symbols indicate numerical values deduced by nu- 
merical differentiation of the results from the Monte Carlo sim- 
ulations [11]. The solid curve represents values calculated from 
the crossover Landau model with u = uo(oo)R~ 4 for R m > 2. 
The dashed curve represents 7^ for R m = 1 corresponding to 
u > 1. 




FIG. 6. Scaled order-parameter <f> = 4>t~^ /Bo as a func- 
tion of — t/t x - The symbols indicate numerical values from 
the Monte Carlo simulations [10,11]. The curve represents 
values calculated from the crossover Landau model with 
u — uo(oo)R~ 4 . 




FIG. 5. Effective susceptibility exponent 7^ as a function 
of t/tx- The symbols indicate values deduced by numerical 
differentiation of the results from the Monte Carlo simulations 
as in Fig. 4. The solid curve represents values calculated from 
the crossover Landau model with u = uo (oo)ir 4 (l + 0.9/r 2 ) 
for 7? m > 2. The dashed curve represents 7^" ff for R m — 1 
corresponding to u > 1. 
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FIG. 7. Scaled susceptibility x~ = X~ T x/^o as a function 
of — t/t x . The symbols indicate numerical values from the 
Monte Carlo simulations [10,11]. The curve represents values 
calculated from the crossover model with u = uo(oo)ii~ 4 . 
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FIG. 8. Effective order-parameter exponent /3 cf j as a func- 
tion of — t/t x . The symbols indicate values deduced by numer- 
ical differentiation of the results from the Monte Carlo simu- 
lations [11]. The solid curve represents values calculated from 
the crossover Landau model with u = iio (oo)R~ 4 for 7i m > 2. 
The dashed curve represents /3 e g for R m — 1 corresponding to 
u > 1. 
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FIG. 9. Effective susceptibility exponent 7^ as a function 
of — t/t x . The symbols indicate values deduced by numerical 
differentiation of the results from the Monte Carlo simulations 
[11]. The solid curve represents values calculated from the 
crossover Landau model with u = uq (oo)R~ 4 for R m > 2. 
The dashed curve represents 7~ ff for R m = 1 corresponding to 
u > 1. 
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